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I. Abstract 

In earlier investigations, the adaptation and implementation of a modified two- level corrections process as 
the onboard targeting algorithm for the Trans-Earth Injection phase of Orion is presented. The objective of 
that targeting algorithm is to generate the times of ignition and magnitudes of the required maneuvers such 
that the desired state at entry interface is achieved. In an actual onboard flight software implementation, 
these times of ignition and maneuvers are relayed onto Flight Control for command and execution. Although 
this process works well when the burn durations or burn arcs are small, this might not be the case during a 
contingency situation when lower thrust engines are employed to perform the maneuvers. Therefore, a new 
version of the modified two-level corrections process is formulated to handle the case of finite burn arcs. This 
paper presents the development and formulation of that finite burn modified two-level corrections process 
which can again be used as an onboard targeting algorithm for the Trans-Earth Injection phase of Orion. 
Additionally, performance results and a comparison between the two methods are presented. The finite burn 
two-level corrector formulation presented here ensures the entry constraints at entry interface are still met 
without violating the available fuel budget, while still accounting for much longer burn times in its design. 

II. Introduction 

U nder nominal operational conditions, the Crew Module’s (CM) 33,361 N main engine performs the 3- 
burn Trans-Earth Injection (TEI) sequence. However, the spacecraft must also have the ability to target 
and execute the maneuvers to return the crew safely to Earth using the backup auxiliary engines. Since the 
total thrust of these eight (8) auxiliary engines is 4,400 N, the burn durations are naturally much longer. 
By way of comparison, that yields a burn duration of 5.5 minutes for the baseline TEI-1 maneuver with the 
main engine, and 55 minutes to perform the same maneuver with the auxiliary engines. 

In the next section of this paper, the formulation of the finite burn algorithm is presented. Next, the 
performance results and a comparison between the two methods is discussed. For all the runs considered, 
the same initial conditions and same entry targets are used. An optimal trajectory will be generated for 
each case to facilitate a comparison with the targeting algorithm results. 

III. Finite Burn Targeting Algorithm 


A. Level I Process 


The Level I process is a differential correction scheme that seeks to identify the transfer arc between two 
position vectors by adjusting the departure velocity. Consider a segment of the trajectory whose endpoints 
(patch points) are denoted by k — 1 and /c, as shown in Figure 1. In a Level I process that employs finite 
burns rather than impulsive maneuvers, the burn arc is considered to be a subsegment of the arc between 
patch points k — 1 and k. The end of the burn occurs at point T. In identifying finite burn arcs, it is 

r i T 


necessary to consider an augmented state vector 


R k V k m k m 9k u k 


, where m k , and u k 
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Figure 1. Level 1 Process 


represent the spacecraft mass, the propellant flow rate, and the thrust direction associated with patch point 
/c, respectively. The goal is to identify a relation between the target, which is the terminal position vector at 
point k ( Rk ), and the control variables. The control variables are the thrust direction ( Uk-i = Uk) and the 
time at the end of the burn (i^)- As in the impulsive formulation, the state transition matrix is partitioned 
into sub-matrices corresponding to each state. The variational equations for each subsegment are then 
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( 2 ) 


for the subsequent coasting subsegment. 

For this formulation, both the initial and final times of the arc (tk-i and tk) are fixed, as are the initial 
position Rk- 1, velocity Vfc-i) and mass rrik- 1- The mass flow rate, m 3fe _ 1 , is also fixed. It is important to 
note that V 7 J~ = Vrf (and therefore SV^ = SVrf), and that 8V^ — = SVf — a^Str + {a? — 

Incorporating these substitutions, the first two variational equations from Equation (1) and Equation (2) 
can be combined to give an expression for 8Rk\ 


SRk = Ak,rGT,k- i$ u k-i + + (a T ~ 


(3) 
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Rearranging this equation into matrix form gives 


r 1 + 

SRk = (A^tGt^-i +Bk,TJT,k-i) Bk,T{a>T ~ ®r) * (4) 

L J OlT 

As in the impulsive Level I method, 1 the minimum norm solution is used to obtain the desired change in 
the control variables such that 

= M T (MM T )~ 1 5R k , (5) 

where 

M = {Akj'Gr.k - 1 + Bk,rJT,k- i) B^t^t ~ ®r) * (®) 

For the first iteration, an initial guess for the duration of the burn is obtained by doing one iteration of an 
impulsive-maneuver Level I correction and plugging the resulting AV value into the rocket equation. The 
initial guess for the direction of the burn is taken directly from the AV vector. 

B. Level II Process 

The Level II formulation for finite burn maneuvers is not conceptually dissimilar from the impulsive maneuver 
scheme. 1 In both cases, the only control variables employed are the positions and times of the patch states. 
As in the impulsive formulation, a minimum norm solution is used to compute the patch point updates. 
However, because the dynamical model for the finite burn algorithm is augmented, the partial derivatives 
are more complex than those in the impulsive maneuver algorithm. The burn is assumed to start with the 
same initial velocity as the terminal velocity of the preceding arc. Thus, the velocity discontinuity occurs 
at patch point /c, the end of the coasting subsegment of the full arc (assuming this arc is then followed by 
another coast arc). Although this problem at first seems identical to the impulsive maneuver targeting, since 
the velocity discontinuity falls between two coast arcs, the partial derivatives for 8V k with respect to 8Rk~ i, 
tk- i, SRki and tk differ due to the thrust segment at the beginning of the arc. 

Extracting 8V k from Equation (2) in the Level 1 formulation section gives 

8V k = Ck,T (SRt — YpStr) + — a^Str) + a k 8tk- (7) 

Recall from the Level I formulation that Vf = V^t at the terminal point of the burn arc and thus that 
8V^ — aj,6tx = SVrf — a^Str + (a^ — a^)8tr- Substituting these expressions into Equation (7) yields 

8V k = Ck,T{8RT — Vt fit t ) + BkjTiifiVrp — a T 8tT) + (a T — a^)StT] 4 - a k 8tk- ( 8 ) 

The first two variational equations embedded within Equation (1) are then summarized as 

(5Rt — VfStx) = A T ,k-i{fiRk-i — y^L^fitk- i) 4- BT,k-i(fiVkLi ~ ^k-ifitk-i) 

+ £ T ,fc-l(<W-l + l) + - fh+^Stk- 1) 

"I" 1 1 1 )> 

(SV^ — a^Str) = Cr,k~i{^Rk~i — VfeLi^fc-i) + — a^^Stk- 1) 

+ + m+^Stk- 1) + / r , fc _i((5m+ - m+ 

+ - ul^Stk- 1). 

For the Level II process, ^ = 0 and u k _ 1 =0. It is still assumed that rh g is a fixed constant, 

i.e. 8m^ k _ 1 = 0. Neglecting these zero terms, and substituting the expressions in Equation (9) into Equation 
(8) , an expression is found for SV k in terms of the state at patch point k — 1 and the state transition matrix, 

SV k = Ck tT [A Tt k-i(8Rk-i - V+^Stk-i) 

4- BT,k-i(fiVkLi ~ ^k-ifitk-i) 4- 4- GT,k-ifiu k _^\ 

4- Bk,T[CT,k-i(fiRk-i ~ V^Stk-i) 4- D T ,k-i{fiVk-i — a>k-ifitk- 1) 

4- HT^k-irii^^fitk-i + Jr.k-ifi^k-i + (®r — Oj,)8tr] + % Stk • 
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In order to write SV r fc - only in terms of the Level II control variables, the first vector equation from Equation 
( 1 ) is used to solve for SV^_^ Su £_ 1} and Stx - Substituting for and SV — a^Stx as shown previously, 
this equation is 

SR k — V k St k = A ki x {SRx — Vx ^Pr) + ^k,r[(SVf — a^Stx) + (a? ~ ( 11 ) 

Dropping the terms that are equal to zero and substituting the same expressions from Equation ( 9 ) that 
were used to form Equation ( 10 ), this first vector equation becomes 

SRk - V k = Ak i T[A.T i k-i{SRk-i — V^Stk-i) + B T ,k - - a^Stk- i) 

+ Erji-irn^ h _^Stk-i + Gx,k-iSu k _^\ 

+ B k} x[Cx,k-i(SR k -i — V^Stk-i) (12) 

+ DT,k-i(SVkL 1 — tik-iStk- i) + H Ti k-i^ k ^Stk-i 

“I - lS^k—i “f" (®r drp^St r p\. 

Further rearranging Equation ( 12 ) into vector form to isolate SV k _ 1 15 Su k _ v and Six yields 


SRk — V k Stk — (A k ,xAx,k - 1 +fifc 1 TC , T,fc-i)(^fc-i — V£i 1 5 tk-i)+ 

(A k ,xBx,k- 1 + Bk,TDT,k-i)at-iStk-i — (A k ,xBT,k - 1 + £fc,T#T l fc-i)ra+_ 1 Jtfc_i 


(A k ,xBx,k-i + B^tDt^-i) (^fc,rGr,jb-i + B^tJt^- i) B^t^t — «t) 


The minimum norm solution is used to solve Equation ( 13 ) for 

T 


Six 


( 13 ) 


Wk-i H-i fa 


in terms of the Level 2 control variables such that 


6u t-i 

Stx 


Here, Z 


= Z T (ZZ T )~ 1 [SR k — V k 6t k 

— {A k ,xAx, k ~i + B k: xCx,k-i)(SRk-i — y^LiStk-i) 

+ (A k ,xBx,k-i + B ki xDx,k-i)^k-iStk-i 

— (A k ,xET,k-l + ^fc,T^T ) fc-l)?hg fe _ 1 ^/t-l]* 


(A k ,xBx,k-i + B ki xDx,k-i) {A k ,xGx,k-i + B k} xJT,k-i) B k ^x{fLx ~ a T ) 


( 14 ) 


. Going back to 


Equation ( 10 ) , this expression is rewritten in terms of the Level 2 controls and the vector ^ SV r j ^_ 1 Su k-1 Stx 
as follows 

SV k = a k Stk + (C k ,TA T ,k - 1 + D k ,TCx,k-i)(SRk-i — V+^Stk-i) 

— (Ck,xB T ,k-i + B kj xBx } k-i)^k-iStk-i + {Ckj'Ex^k-i + Dk^Hx^-i^g^Stk-i 


+ [ (C k ,TB T ,k - 1 + Dk,xD T ,k-i) (C k ,TGx,k - 1 + L)fc 5 T«/T,fc-i) B> k ,x{a T ~ a x) 

r n T 1 

Substituting the expression for 


H-i 

Stx 


( 15 ) 


SV^ Su ^_ 1 Stx J into Equation ( 15 ) gives 

SV k = a k St k + (G k ,TAx,k - 1 + B k ,xCx,k-i){SRk-i — V^Stk-i) 

— (Ck,TBx,k- 1 + Dk,xDx,k-i)a k -iStk-i + {G k ,TEx,k-i + Dk^Hx^k-i^g^Stk-i 
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+ D^tDt^-i) (C^tGt^-i + B^tJt^-i) D^t^t ~ ®r) Z T (ZZ T ) l . (17) 
Then, rearranging Equation (16) to isolate the control variables yields 

SV^~ = [(Ck,TA T ,k - 1 + Dk t TCr t k- 1) _ Z{Ak^TATfr - i + ^.tCt^-i)]^-! 

+ ([(Cfc,T^T,fc-l + A;,T#T,fc-l) - Z(Ak,TE T ,k - 1 + ^fc,T^T,fc-l)]m+ 

— [(Cfc } T#T 5 fc-l + Afe.TAr.fc-l) — Z(Ak,rBT,k - 1 + ^fc,T^>T,A;-l)]a^_ 1 (18) 

— [(C^t^t^-i + T^fc^C'^fc-i) — Z{Ak,TA T ,k-i + Sfc,rCt,fc_i)]V]^l 1 )Jtjb_i 

+ ZSRk + (cifc — )^fc* 

Thus, the partials of are 


d5Rk-i 


= (C^tAt^-i + D^tCt^-i) — ^(A^tAt^-i + ^tCt^-i), 


: = [(^fc,T^T,fc-l + Ac 5 T#T,fc-l) — Z(Ak,TET,k - 1 + ^ 

— [(Cfc,T#T,fc-l + ^fc 5 T^T,fc-l) “ (Ak t T B t , k—l + ^fe } T^>T,fc-l)]«fc_i 

— [(Cfc,T^T,fc-l +^fc 5 TCT,fc-l) — ^(^fc,T^T,fc-l + ^fc J TC'T J fc-l)]Vfcti, 


(% -^”)* 


Because it is assumed that the arc from patch points fcto/c+lisa coast arc, the partial derivatives of SV^~ 
are the same as for the impulsive case. The partial derivatives of AV*. are therefore 

,, 0 — — = —(C^tAt^-i + B^tCt^-i) ~ %{Ak,TA T ,k - 1 + i)> 
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The desired changes to the control variables are found using the minimum norm solution, 
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C. Maneuver Sum Constraint 


In addition to the velocity continuity constraint, several trajectory constraints are imposed in the Level II 
process. 2 One such constraint is on the total AV sum of the maneuvers. The finite burn formulation of this 
constraint is based on the impulsive maneuver sum constraint, 3 and in fact uses the exact same approach to 
place the relevant partial derivatives into the M matrix. Only the composition of those partial derivatives 
and the error calculation process change. 

Denote the patch point where the burn occurs as patch point /c, rather than k — 1 as in the previous 
sections. This is so that the Level II control variables for this constraint are still R k -u tk-u Rk , Rk+ i } 
and t/c+i, as is the convention for this process. To derive the burn maneuver constraint, it is necessary to 
determine the partial derivatives of the magnitude of A V*., the AV resulting from the burn maneuver at 
patch point fc, with respect to the Level II control variables. Prom the rocket equation, AVfc is given by 

= -I, p90 ln( 1 - ThgfcAftarn ), (22) 

m k 

where A tbum = tr —t k - The partial derivative of AVk with respect to the burn duration at patch point k 
is given by 

)(^ 

(J Alburn Album ?Tlk 

Next, the partial derivatives of A tbum with respect to the control variables are necessary. Prom the definition 
of A tbum’, h is evident that an expression for Str in terms of the control variables is required. This is 
determined using the variational equations from points k — 1 to fc, k to T (the termination of the burn 
segment), and k + 1 to T; 
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(24) 


(25) 


(26) 


Recall that V? = V T , and therefore SV^ = SV T . Then the second vector equation of Equation (25) is 
subtracted from the second vector equation of Equation (26) to isolate Str 


(■ O r - a+)St T = -C T ,k( SR k - V^Stk) - D Tik {6V+ ~ a%5t k ) 

— H T , k {8m k + ihg k &tk) “ Jr,kSu k + CT,k+i{SRk+i ~ ( 27 ) 

+ D r,k+ 1 (^+ 1 . - a^ +l St k+ 1 ). 

Since there is no velocity discontinuity at patch point fc, SV^ = SV^\ thus, from Equation ( 24 ), 

<^fc + - St k = {W k ~ - d k St k ) + (d k - d k ) 

= C k ,k-i(SRk-i ~ V+-i St k-i) + R>k,k- i^V^i - d^dtk- i) + (a k - a^)5t k . 


Substituting this expression back into Equation (27) yields 

(«t — a^)Str = —CT, k (SR k ~ Y k ^k) — D T ,k [Ck,k-i(5Rk-i — V^Stk-i) 

"h ^k— l^k— l) T ( & k &j^)] HT,k(&1Tlk 

— JT,k^t + Cr,fc+i(&Rfc+i — V k ~ +1 5t k +i) 

+ D T ,k+ l^fc+i ” Afc+i^fc+i)- 
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The term SV k ~_ 1 — a^^Stk- 1 may be written in terms of the control variables by rearranging the first 
variational equation of Equation (24), 

SV+_1 — = B kk _ i[(SB k — V k St k ) — — V k -\St k )\. (30) 

Thus, Equation (29) may also be written as 

(' a t ~ 4)^t = — Cr 5 fc(^fc — ^ + <^fc) — D T ,k\Ck,k-i{$Rk-i — 

+ — n ^fc) — Ak,k - 1 (SRk-i — Vk—i St k )] + — 

— H T ,k(firrik + St k ) — Jr,kSu k + CT 5 fc+i(^fc+i — V k+1 5t k +i) 

+ D Tik+1 (SV k - +1 - a k+1 St k+ 1). 

Grouping like terms and rearranging to combine the non-control variables into a single vector yields 


(flji Cip)St r p 


= — (Cr,fc + Ar,fcAfe,fc-i^4i)(^fc — Vk^tk) 

- D T}k (Ck,k-i ~ Dktk-iBktk-iAktk-^ifiRk-i - V^Stk-i) 
~ D T,k(H ~ n k) St k ~ H T ,k( Sm k + m^Stk) 

+ Ct,fc+i(&Rfc+i — Vfc+i^fc+i) ~~ Dp ik +ia k+1 St k +i 


+ “ Jr,fc At,/c+1 


It has been shown, in the derivation of impulsive targeting algorithm, that (C k ,k-i — D kik - iB k k-i A k , k -i) 
fc ; thus, Equation (32) can be somewhat reduced to 

(a^ — 4 )^t = — (Cr,fc + ^fc) 

- ^T 5 /c(d^ - “ H T , k (6m k + m+ <ft fc ) 

+ Ct,fc+l(£Rfc+l — Vfc+i^fc+l) — ^T,fc+l^fc + l^fc+ 1 


+ — «^r,fe ^T,fc+i 


r _ -| t 

In order to express 8u k SV k+1 in f erms of the control variables, the first vector equations from 

Equation (25) and Equation (26) are used. Because Vp = Vp , it is also true that 5Rt — Vp = SRt — 5Vf. 
Thus, 

A Tik (8R k ~ V r k 8t k ) + Bp jk (8V^ - a^St k ) + E T}k (Sm k + rh+ k 5t k ) + G?T,fc^4 
= 4r,fc+i(£Rfc+i - V^+i^fc+i) + -®r,fc+i(^+i “ Ofe+i^fc+i)- 
Using the relationship in Equation (30), the previous equation may be rewritten as 

AT,k($Rk ~ V k Stk) + B T ,k[Ck, k -i(SR k -i — Y k _iSt k ~\) 

+ ^ k , k -iB kk _ 1 \{5R k — V k 5t k ) — A k ^ k -i(5R k -i — V k _^ St k )\ + ( a k — a£)5t k ] 

-4- (35) 

+ B T , k ( Sm k + rhg k St k ) + G Ti k$u k 

= A Tik +i(5R k +i - V k+1 5t k+ 1 ) + B T ^ k+1 {5V kJrl - a k+1 St k+1 ). 

Rearranging this equation yields 

{Ar, k + B T ^ k D k ^iB k 1 k _ 1 )(5R k — V k St k ) + Bp ik (C kjk - i 

- B kik _ 1 B k l_ 1 A k ^i)(5R k -i - V+^Stk) 

+ BtM^ ~ a%)8t k + E T}k (5m k + rh+ k 6t k ) 

- _ (36) 

— ^.T,fc+i(<^fc+i — V k+1 5t k + 1 ) + B T > k+1 a k+1 5t k +i 
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Again using the relationship (Ck,k - 1 — k-iAk,k-i) = B k _ x k , Equation (36) becomes 

(Ar,k + B T ,kDk,k-iB k ^ k _- L ){5Rk - V k 5tk) + BT,kB k ]i lk (6Rk- 1 - V^Stk) 

+ B T ,k(a k - a k )5tk + E T ,k(5mk + rhg k 6tk) 

- AT,k+i(dRk+i ~ + B T ,k+ia k+1 6tk+i (37) 


r 

Sut 

—Gx,k B Ti k+ i 

k 

. 6V k + i . 


Letting S = 

[ SV k+1 


—Gr,k &T,k+ 1 

iT 

in terms of the control variables: 


, the minimum norm solution to Equation (37) gives an expression for 


™k + i 


= S T (SR r )~ 1 [(A T>k + B^kDk^B-l^iSRk - V k 6t k ) 

+ BT,kB k \ k (5Rk- 1 — VfeLi^tfe) +-Bj' ) fc(a fc — a k )5tk 
+ B T ,k(Smk + rhg k 6tk ) — A Tl k+i(SRk+i — V^~ +1 5tk+i) 

+ B T ,k+ia k+1 6tk+i]- 


(38) 


This can be substituted back into Equation (33) to obtain an expression for Str in terms of the control 


variables. Let S ■ 


~Jr,k Ar,fc+i 


S T (SS T ) 1 ; Equation (33) can then be written as 


(a T — a^)Str = — (Cr,fc + ^T,fcAfe,fc- i^k,k-i)(^^ k ~ Vk^k) 

— D Tik B k _ lik (5R k .i — — DT^ki&k ~ 

— H Tf k(^k + fh^Stk) + C f T,fc+l(^fc+ 1 - Vfc+i^fc+l) 
— &T,k+ lQ>k+l$tk+l 

+ S[(A T ,k + B T ,kBk,k-iB k 1 k _ 1 )(SRk — V k ~ St k ) 

+ ^T i kB k _ lk (SRk-i — V^Stk) + BT,k(^k ~ 

+ ET,k(fi m k + "hJk^fc) ” A.Tfr+l{SRk+l — V k+1 6tk+l) 

+ #T,fc+ lOfc+i^fc+l]- 


It should be noted that the term <5mfc is still present in this expression; this term is a function of the control 
variables at previous maneuver patch points and will be dealt with at the end of this derivation. Combining 
like terms reduces Equation (39) to 

(a T — a^)StT = — (Dr,kB k \ k — SBT,kB k \ k )(5Rk-i — 

— [{Cr,k + Dr^Dk’k-iB^l^) — S(A T ,k + — V k Stk ) 

— (B>r,k ~ SB T ,k)(^k ~ )^fc — (^T,k — SE T ,k)(Srrik + ra+Jt k ) (40) 

+ (Cr 5 fc+i — SA T ,k+i)(SRk+i — V k+1 8tk+\) 

— {D Ti k + 1 — &B T ,k+i)ak+i^ k + lm 

Finally, let A olt = — a^. This term can be written in terms of its magnitude and direction; 


A aT = 1 1 Aa^l | 


A aT 
||Aot|| 


(41) 


At point T, Aar is equal to the applied thrust. 


Thus, the unit vector ||^^|| is equal to the unit thrust 


8 of 16 


American Institute of Aeronautics and Astronautics 



direction u k = tt^tt • Equation (40) can therefore be written as 

ll^fc II 

\\Aa T \\ukSt T = -(D Tik B^ i,k ~ SB Ttk Bj-\ ik )( 5 Rk- 1 - St k ) 

- [(Cr,fc + - S(A T ,k + St k ) 

- (Dr,k ~ SB T ,k)(a k ~ Q> k )St k - (H Tik - SE Tjk )( 5 m k + ra+ St k ) 

+ (Cr,fc+ 1 - SA T ,k+i){SRk+i — V k+ i^fc+i) 

- (Btj c+i — SB T , k -\-i)a^ +1 St k+ i. 

Dividing both sides of Equation (42) by ||Aar||, 

tutor = IIAorll - SB^B^ ^R^ - V+^SU) 

— [{Ct , k + ^T,fc^)fc,fc-i^fc_i) — S(Ar, k + ^T,fc^fc,fc-i^fc_i)](^fc — V k St k ) 

— (B>r,k — SB Tjk )(a^ — a^)St k — {H T ,k ~ SE T:k )(Sm k + m+ St k ) 

+ ( Cr,k+i — SAT )k +i)(SRk+i — V k +iSt k +i) 

— (Ar,fc+ 1 — SBT.k+l^k+ltok+l]- 

Because u k is a unit vector, both sides of Equation (43) can be pre- multiplied by u £ to give 

tor = ]^\[-( D T^B k \ k - SB T j,Bj£ l 9 k )( 8 R k - 1 - V+_ , 6 %) 

— [{pT,k + Br^D^k-iB^l-i) — S(A T ,k + ^T,fc^fc 5 fc-i^fc_i)](^fc — Vjj~ St k ) 

— ( B T ,k — SB T ,k){a k — & k )Stk — (Hr,k ~ SE T ,k){Sm k + 

+ (Cr,fc + i — S^4t,/c+i)(^/c+i — Vfc+i^fc+i) 

— (-^Tjfc+i — SBT,k+i)& k +iSt k +i]. 


Using Equation (44) and recalling that A tbum n 
the control variables are 


tr—tk, the partial derivatives of A tbum n with respect to 


dAtburn n 

dRk-i 


||Aot|| 


(D T ,kB k \ k - SBr,kB k K k ), 


~i£T = pfcn (DT '‘ B ‘-‘-’‘ ~ 

9A tburn n _ U k . d-1 \ at 


dR k 

dAtbum n 

dt k 


||Aa T || 


\{pT,k + D T ,kDk,k-iB k fc _ 1 ) — S(A T>k + B T ,kDk,k-iB k 


i([(Cr,fc + D T ,kDk,k-\B k I_ x ) - S(A T ,k + B T ,kDk,k-iB k l_ 1 )]V k 


dt k ||Aa T |r -*,«-«,K~L-k,k- is 

— {Dr,k — ^B T ,k){a k — a k ) — (H T ,k ~ SET,k)'dig) — 1, 

dAtburn n _ U k . , r, A \ 

~dR^r - iiiMi (CT ’‘ +i “ T ' i+l) ’ 

SA tburn n _ ^k r (r<_. . O A - I ( TA . . CD_ 


||Aa T || 


[(Cr,k + 1 — &Atjz+ i)V k +i + (Br,k+i — SB Tjk+ i)a k 


The term containing Sm k can be expressed as a function of the change in the initial mass, Jmo, and the 
burn durations of all maneuvers prior to patch point fc, 


Sm k — StTIq ^ ^ Wig ^tburrij j 
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where j denotes the maneuver number and n is the maneuver number at patch point k. Because the initial 
mass is fixed, this equation reduces to 

71 — 1 

5m k = ^ ^ TfLgAt} )urn ^ . (47) 

3 = 1 

Since m g is a fixed, constant value, Srrik is a function only of the previous burn durations. Thus, A t burrin has 
a dependence on the positions and times associated with any previous maneuvers that have occurred. Also, 
AVfc has a direct dependence on the mass at the beginning of the burn (and therefore on the previous burn 
durations), so that must be taken into account as well. Using the chain rule, the final form of the partial 
derivative of the constraint a (the sum of all the burn A Vs) with respect to any control variable f3 k in the 
set of control variables associated with patch point k is 


da _ dAV n dAt burrin dAVj dA t burrin 

d(3 k dAt burrin d(3 k dAtb urn n dp k 


where 

dAVj _ dAVj dA Vj dAt hurnj 

dAtburn n dAtb urn n & ^^burrij dAt burn n 


(49) 


These partial derivatives are then placed into the M matrix in accordance with the existing maneuver 
constraint formulation. 


IV. Simulation and Results 

The two level targeting algorithm is applied in the following section in both the impulsive and finite burn 
configurations. This is done in order to demonstrate its performance under both configurations and to gain 
insight as to applications of each. Then, as a final metric, an optimal trajectory is generated to use for 
performance comparisons. 

For each case, the same initial conditions of the lunar orbit will be used as follows: 

• Epoch: 4- Apr-2024 15:30:00 TDT 

• Initial mass: 20339.9 kg (total fuel = 8063.65 kg) 

• Main Engine Thrust: 33,361.6621 N 

• Main Engine Isp: 326 sec 

• Auxiliary Engine Thrust: 4,448.0 N 

• Auxiliary Engine Isp: 309 sec 

• State (J2000 Moon-centered inertial frame): 

- X: -1236.7970783385588 km 

- Y: 1268.1142350088496 km 

- Z: 468.38317094160635 km 

- Vx: 0.0329108058365355 km/sec 

- Vy: 0.589269803607714 km/sec 

- Vz -1.528058717568413 km/sec 

Likewise, the same terminal target conditions will be used for each case as shown below: 

• Geodetic Altitude (km) 121.92 

• Longitude (deg) 175.6365 

• Geocentric Azimuth (deg) 49.3291 

• Geocentric Flight Path Angle (deg) -5.86 


10 of 16 


American Institute of Aeronautics and Astronautics 



A. Impulsive Algorithm Example 

The first case is representative of a nominal Earth return during which the maneuvers are performed by the 
CM main engine. For this case, the impulsive targeter is executed and the initial guess file consists of 12 
” patch points” or states taken from an optimized trajectory. The first patch point corresponds to a state ad 
epoch 20 minutes prior to the optimized TEI -1 TIG. The interior patch points correspond to the states and 
epochs at each of maneuver locations (TEI-1, 2, 3 and TCM 1, 2, 3) and some additional waypoints along 
the trajectory. The final patch point in the initial guess is the state and epoch at the desired entry interface 

(EI). 

For this case the impulsive targeter is executed in order to find a feasible trajectory that satisfies the 
specified terminal constraints, while minimizing the AV sum of the individual maneuvers. The targeter 
converges on a feasible solution after 15 iterations and 2.34 minutes of computation time (Figure 2). The 
total AV of this solution is 1.3484 km/s. The AV for each maneuver is given in Table 1. 
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Figure 2. Initial Guess and Impulsive Solution 


Table 1. Maneuver Data 


Maneuver 

Ay (km/s) 

TEI-1 

TEI-2 

TEI-3 

0.6753 

0.2617 

0.4114 


B. Finite Burn Algorithm Example 

Assuming the use of the main engines for all three maneuvers, the finite burn algorithm is able to converge 
on a final solution in 24 iterations, with a total AV of 1.3353 km/s. The maneuver sum constraint is active 
for this case, with a maximum allowable Ay of 1.40 km/s. The full solution is shown in Figure 3. Figures 4, 
5, and 6 show closeup views of TEI-1, TEI-2, and TEI-3, respectively. The solid red sections in each of these 
closeups indicate the thruster burn segments of each arc. The burn duration, propellant mass consumption, 
and AV for each maneuver are given in Table 2. The total computation time required to achieve a solution 
is 6.57 minutes. 
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Finite Bum Trajectory (Main Engines) 
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Figure 3. Initial Guess and Finite Burn Solution with Main Engines, Full Trajectory 
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Figure 4. Initial Guess and Finite Burn Solution with Main Engines, TEI-1 Closeup 


Table 2. Burn Data Using Main Engines 


Maneuver 

Duration (s) 

Prop. Mass Consumed (kg) 

AV (km/s) 

TEI-1 

377.6462 

3939.530 

0.6546 

TEI-2 

131.2084 

1368.740 

0.2650 

TEI-3 

184.1434 

1920.947 

0.4157 


12 of 16 


American Institute of Aeronautics and Astronautics 


Finite Bum Trajectory (Main Engines) 


Figure 5. 


Figure 6. 
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Initial Guess and Finite Burn Solution with Main Engines, TEI-2 Closeup 
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C. Finite Burn Algorithm with Auxiliary Engines 

For this case, the maneuver sum constraint is relaxed to a maximum AV of 1.50 km/s. Using only the 
auxiliary engines, the algorithm converges in 18 iterations, with a total AV of 1.500 km/s. As in the previous 
section, the solid red portions of the final trajectory closeups indicate the segments of the trajectory in which 
the engines are thrusting. Figure 7 shows the impulsive initial guess and the final finite burn solution for the 
full return trajectory. Figures 8 through 10 show the closeup views of TEI-1, TEI-2, and TEI-3. Burn data 
for each maneuver is listed in Table 3. The solution is identified after 7.67 minutes of computation time. 
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Figure 7. Initial Guess and Finite Burn Solution with Auxiliary Engines, Full Trajectory 


Table 3. Burn Data Using Auxiliary Engines 


Maneuver 

Duration (s) 

Prop. Mass Consumed (kg) 

AV (km/s) 

TEI-1 

3220.3465 

4725.536 

0.8039 

TEI-2 

934.8663 

1371.823 

0.2796 

TEI-3 

1242.5484 

1823.316 

0.4165 


D. Optimized Finite Burn Trajectory 

As a final step in this analysis, the previous trajectory generated from the finite burn targeter with auxiliary 
thrusters was optimized using Copernicus. Again, the same initial conditions and initial guess file were used 
to initiate both runs. The results of this optimized run are available in Table 4. The total AV for the 
optimal run is 1.2413 km/s. This is an improvement of approximately 0.25 km/s of AV over the finite burn 
targeting solution. 


Table 4. Optimal Burn Data Using Auxiliary Engines 


Maneuver 

Duration (s) 

Prop. Mass Consumed (kg) 

AV (km/s) 

TEI-1 

3327.008 

3501.37 

0.6040 

TEI-2 

873.040 

1504.85 

0.2993 

TEI-3 

1312.285 

1828.22 

0.4059 
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Figure 8. Initial Guess and Finite Burn Solution with Auxiliary Engines, TEI-1 Closeup 
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Figure 9. Initial Guess and Finite Burn Solution with Auxiliary Engines, TEI-2 Closeup 
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Finite Bum Trajectory (Auxiliary Engines) 
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Figure 10. Initial Guess and Finite Burn Solution with Auxiliary Engines, TEI-3 Closeup 


V. Conclusions 

This paper presents a two-level targeting algorithm for finite burn maneuvers. The algorithm is a modifi- 
cation of the impulsive targeting scheme described in Reference 2. Formulations for the finite burn versions of 
both the Level I and Level II process are discussed, and a total mission cost constraint, originally developed 
for the impulsive algorithm, is modified and adapted to the finite burn problem. The algorithm is tested for 
the trans-Earth injection phase of the Orion mission. Results are compared for three different cases: using 
the impulsive algorithm, using the finite burn algorithm with main engines for all three maneuvers, and using 
the finite burn algorithm with only the auxiliary engines for the final two maneuvers. Results show that 
the finite burn algorithm is able to converge on a feasible solution even in the case of a main engine failure 
following TEI-1. Due to the additional complexity of the finite burn equations, the algorithm requires more 
computational time than the impulsive targeter, but it still shows significant improvement over an optimizer 
in terms of the computation cost. 
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